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ABSTRACT 

A  new  theoretical  approach  to  the  problem  of  radar- 
target  scintillation  was  developed.  It  was  applied  to  slender, 
axially  symmetrical  targets  that  are  much  longer  than  the  RF 
wave  length.  The  radar  scintillation  is  analyzed  and  described 
in  terms  of  statistical  parameters.  The  following  sets  of  quantities 
were  derived  and  computed: 

(a)  mean  radar  cross  sections,  their  RMS  fluctuations  and 
average  lobing  frequencies; 

(b)  mean  target  centroids,  their  RMS  deviations  and  average 
meandering  rates. 

These  and  other  variables  were  expressed  as  functions 
of  the  aspect  angle,  radio  frequency  and  of  configurational 
detai Is. 
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I.  INTRODUCTION  TO  THE  PROBLEM  OF  RAD 4 R -TARGET  SCINTILLATION 

The  radar  characteristic  known  as  "target  scintillation" 

was  defined  by  Muchmore  [11  as  a  physical  phenomenon  that  consists 

of  two  effects:  (1)  the  fluctuation  of  the  rectified  amplitude  of 

the  target  echo  as  observed  in  the  video  circuits,  and  (2)  the 

wander  of  the  apparent  centroid  of  target  position  as  generated 

* 

by  the  output  of  a  tracking  circuit.  The  scintillation  of  the 
stars  also  consists  of  two  effects:  fluctuation  of  magnitude  and 
variation  of  position.  However,  it  was  pointed  out  by  Muchmore 
that  the  causes  of  stellar  and  radar  target  scintillations  are 
distinct.  Star  light  twinkle  is  caused  by  the  disturbances  of 
the  intervening  medium  (atmosphere)  whereas  radar  aircraft 
fluctuation  is  produced  by  the  interference  between  many  scattering 
points  of  a  usually  complex  target.  It  will  be  shown  later  that  the 
properties  of  a  radar  receiver  also  may  affect  the  character  of 
radar  scintillation. 

For  the  design  of  an  efficient  radar  and  weapons  system, 
it  is  important  to  obtain  a  thorough  understanding  of  the  analytical 


In  the  present  report  we  are  concerned  with  target -centroid  wander 
in  the  angular  space  of  the  radar  field  of  view.  However,  similar 
methods  may  also  be  applied  to  circuits  producing  radar  range  and 
radial  range  rate  (Doppler).  We  shall  follow  Muchmore's  definition 
of  "target  scintillation"  although  some  authors  deviate  from  this 
nomenclature.  For  instance  Barton  (Ref.  [27])  refers  to  "amplitude 
noise"  as  scintillation  which  other  authors  call  "fading  noise." 

The  term  "glint"  is  generally  applied  to  "angular  noise,"  the 
quantity  "radar  jitter"  is  a  collective  term  that  also  includes 
the  effects  of  receiver  and  servo  noise. 


** 


See  Appendix,  Secs.  A. IV  and  A.V.3. 
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and  statistical  character  of  the  target  echo.  The  latter  enters 
into  two  phases  of  radar  operation.  During  the  stage  of  target 
detection  and  acquisition,  the  mean  radar  cross  section  (abbreviated 
"RCS"  in  the  sequel),  its  fluctuation  and  time  rate  of  fading  play 
an  important  role.  The  success  of  this  operation  can  be  measured 
and  computed  by  means  of  tables  and  graphs  of  pertinent  references 
(cf.  [2],  [3],  [4]).  An  important  entry  of  these  tables  is  the 
"category  of  the  target,"  which  is  specified  by  several  parameters. 

One  important  characteristic  is  the  average  duration  of  the  fading 
interval  as  compared  to  the  pulse  repetition  and  scanning  periods 
of  the  radar. 

During  the  target  tracking  and  missile  semi -active 
homing  stages  of  a  missile -target  intercept,  similar  problems 
relative  to  radar  performance  are  involved.  The  RMS  fluctuation 
of  the  mean  target  centroid  and  its  rate  of  oscillatory  meander 
are  both  significant  and  somewhat  parallel  in  concept.  The  problem 
of  "glint"  caused  by  the  erratic  wandering  of  the  centroid,  especially 
for  extended  targets  at  close  ranges,  is  akin  to  the  multiple- 
target  problem  which  has  not  yet  been  solved  in  a  satisfactory 
manner. 

The  knowledge  of  hostile  target  scintillations  can  be 
obtained  (a)  by  actual  flight  observations,  (b)  by  static  radar- 
ground  tests  carried  out  on  models  and  (c)  by  analysis  and 
computation,  provided  the  configuration  of  the  enemy  craft  is 
known  from  photos  or  other  means  of  intelligence.  The  present 
analysis  belongs  to  the  latter  category.  It  offers  the  following 
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advantages:  (1)  cost  reduction  by  avoiding  hardware  and  test 

operations  during  preliminary  design  stages  of  the  defensive 
weapon  system,  (2)  flexibility  in  the  choice  of  systems 
parameters,  and  (3)  mathematical  evaluation  of  possible  schemes 
for  dealing  with  heavily  fluctuating  targets.  These  schemes 
include  the  selection  of  polarization,  radio  frequency  modu¬ 
lation,  narrow-band  filters  the  center  frequencies  of  which 
can  be  adapted  to  measured  instantaneous  frequencies,  adjustable 
automatic  gain  controls,  etc. 


-  3  - 
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II.  BACKGROUND  OF  PAST  AND  PRESENT  RESEARCH 
ON  THE  RADAR  SCINTILLATION  PROBLEM 

Various  distinct  approaches  to  the  problem  of  radar 
target  scintillation  can  be  found  in  the  literature.  The  experts 
in  the  field  of  radar  signatures  (K.  M.  Siegel  and  collaborators, 
see  Ref.  [5]),  developed  an  extensive  physical  and  mathematical 
theory  by  means  of  which  the  RCS  of  a  complex  structure  (consisting 
of  hundreds  of  stations)  can  be  computed  as  a  function  of  aspect 
angle,  frequency  and  polarization.  Good  agreement  between  analytical 
results  and  measured  radar  echoes  is  claimed.  Some  of  the  under¬ 
lying  assumptions  and  mathematical  principles  will  be  discussed 
in  Appendix  A. I.  The  statistical  results  of  this  research  team 
are  limited  to  concepts  of  so-called  "first-generation  statistics," 
such  as  a  (mean  cross  section) ,  S  (standard  deviation  of  a)  etc. 

No  use  is  made  of  quantities  appearing  in  "higher-dynasty  statistics" 
such  as  multiple  joint  distributions,  correlation  functions  ard 
spectra. 

A  different  type  of  a  statistical  theory  of  aircraft 
scintillation  was  developed  by  Delaro  [6]  and  Muchmore  [1], 

Their  results  appear  in  terms  of  spectra  and  probability  densities. 
However,  their  theory  applies  to  simple  airframe  structures  that 
can  be  considered  as  an  assembly  of  a  few  individual  scatterers 
or  of  a  few  simple  areas  of  uniformly  distributed  scatterers. 

Related  problems  were  analyzed  in  Refs.  [7],  [20],  [211,  and  [22], 
Some  interesting,  statistical  work  was  carried  out  by  Robert  W. 
Kennedy,  Major  USAF.  He  was  primarily  concerned  with  radar 
cross  sections  of  simple,  satellite -type  targecs.  (See  Ref.  [8]). 

A  natural  extension  of  the  above-cited  work  would  call 
for  a  combination  of  the  ideas  incorporated  separately  in  the 
investigations  of  the  aforementioned  research  teams.  It  would 
then  yield,  for  instance,  the  RCS  spectra  of  complicated  target 
configurations.  The  preparation  of  such  quantities  presents  a 

-  4  - 
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cumbersome  computational  chore  but  it  can  be  accomplished. 

However,  the  results  would  still  be  inadequate  for  the  design  of 
radars  against  specific  complex  targets,  for  the  following  two 
reasons : 

(1)  "Spectra,”  in  a  strict  sense,  are  only  applicable  to 
stationary  processes.  They  furnish  the  mean  square  amplitude 
densities  in  the  frequency  domain,  averaged  over  an  infinite 
interval  of  time.  The  radar  operators  in  the  field,  however,  are 
well  aware  of  the  fact  that  the  signatures  of  slender  targets 
vary  markedly  in  magnitude  and  spectral  content  when  turning  from 
broadside  to  frontal  aspect.  What  is  needed  here,  is  the  knowledge 
of  temporary  signal  behavior.  The  concept  and  the  analysis  of 
"locally  stationary  processes"  will  be  introduced  and  utilized 

in  the  present  report  to  cope  with  this  phase  of  the  problem. 

(See  Appendix,  Sec.  A. II). 

(2)  Often  one  deals  with  complex  targets  consisting  of  many 

component  scatterers  that  are  somewhat  randomly  distributed  over 

the  length  of  the  body.  Radar  interference  between  these  scatterers 

* 

generates  many  nulls  which  appear  randomly  distributed  in  time. 

In  the  extreme  case  of  a  Poisson  distribution,  these  lobes  generate 
a  flat  spectrum  according  to  Ref.  [9].  This  very  fact  lies  at  the 
heart  of  the  multiple-target  problem.  If  one  deals  with  narrow- 
band  targets  in  the  background  of  wide-band  noise,  one  can  design 
narrow-band  filters  that  emphasize  the  target  signals  and  wash  out 
the  noise.  If  one  deals  with  wide-band  target  signals,  no 
ordinary  filter  is  helpful.  If  one  washes  out  the  noise  one 
filters  the  target  signal  as  well.  The  physical  implication  of 
this  fact  is  well  known  to  the  radar  operator  in  the  field.  If 
he  tracks  two  targets  of  equal  strength  appearing  simultaneously 
and  a  fraction  of  a  beam  width  apart,  the  radar  axis  will  wander 
erratically  from  one  target  to  the  other  and  sometimes  stray 

* 

even  if  the  target  turns  at  constant  angular  rate. 

In  many  references  "target  signals"  are  separated  into  two  parts: 

(1)  "true  target  signals"  such  as  "true  angular  positions"  and 

(2)  "target  noise"  such  as  "glint."  For  other  quantities  such  as 
the  "RCS"  a,  this  separation  is  not  as  simple.  For  some  variables 
of  closed-loop  missile  guidance,  +his  separation  becomes  even 
arbitrary  and  depending  on  the  ch. ice  of  the  analyst. 

-  5  - 
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beyond  these  points  until  the  signal  gets  lost.  If  the  receiver 
applies  a  narrow  band  tracking  filter  or  a  slower  antenna  servo, 
the  radar  axis  still  undergoes  excursions  of  the  same  magnitude, 
only  at  a  slower  rate.  Apparently  what  is  needed  here,  in  order 
to  cope  with  this  difficulty,  is  a  better  knowledge  of  the  so-called 
"instantaneous  frequency"  and  of  other  parameters  developed  by 
S.  0.  Rice  under  the  heading  of  "Statistics  of  the  Count  of  Zero- 
crossings.”  (Ref.  [10]).  J.  W.  Follin,  Jr.,  of  APL  proposed  to 
apply  this  segment  of  Rice’s  random  noise  analysis  to  the  problem 
of  radar-target  scintillation.  The  present  author  performed  the 
mathematical  analysis  connected  with  this  problem  and  the  Conductron 
Corporation  carried  out  the  numerical  computations,  (cf.  References 
[11],  [12],  [13],  [15],  [19]). 

The  major  results  of  this  investigation  will  appear  in 
Section  VII.  For  our  specific  task  a  particular  type  of  target 
configuration  was  chosen  and  used  in  our  numerical  calculations. 
Tactical  geometry  and  the  choice  of  target  coordinates  will  be 
explained  in  Section  III.  General  definitions  and  symbols  are 
given  in  Section  IV,  specific  definitions  and  refined  statistical 
notation  in  Section  V.  The  underlying  assumptions  which  primarily 
pertain  to  our  special  project,  are  compiled  in  Section  VI.  In 
addition,  we  used  a  number  of  general  mathematical  principles  which 
were  very  helpful  in  rendering  the  computational  chores  tractable. 

In  order  not  to  encumber  the  main  stream  cf  information  contained 
in  the  major  body  of  this  report,  mathematical  principles  and 
pertinent  theorems  were  relegated  to  special  sections  of  thr 
appendix,  as  were  refined  details  of  the  mathematical  formulation 
of  the  problem  and  of  the  derivation  of  our  formulae. 
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III.  TACTICAL  GEOMETRY  AND  TARGET  COORDINATES 

The  tactical  geometry  and  target  coordinates  are  explained 
in  Figure  1.  Two  cases  must  be  considered: 

(1)  Monostatic  Case 

This  case  refers  to  the  operation  of  a  tracking  radar 
and  the  action  of  an  actively  homing  missile.  The  radar-line  of 
sight  and  the  target -longitudinal  axis  form  a  plane  (P)  which  is 
depicted  in  Figure  1.  Of  special  interest  are  vehicles  with 
axial  symmetry.  For  this  special  configuration  the  RCS  problem 
becomes  two-dimensional  since  one  spatial  dimension  (Z-coordinate) 
can  be  suppressed. 

The  choice  of  the  coordinate  system  is  shown  in 
Figure  1.  The  origin  coincides  with  the  nose  tip  (station  "0") , 
the  x-axis  with  the  longitudinal  axis  of  the  target,  and  the 
y-axis  lies  in  plane  P.  The  vehicles  are  divided  into  N  sections 
whose  end  points  on  the  contour  or  some  other  convenient  locations 
are  selected  as  the  stations  which  represent  the  individual 

component  scatterers. 

More  details  which  are  important  in  the  analysis  are 
sketched  in  Figure  2. 

Besides  the  stations  ,  the  positions  of  the 
instantaneous  centroid  S  and  of  the  mean  centroid  5  are 
important.  In  References  Tl]  and  [6],  the  latter  are  designated 
as  the  "apparent"  and  "effective"  radar  centers  respectively. 

It  is  convenient  to  refer  the  so-called  error  angles  to  the 

line  of  sight  toward  S.  Then  one  has  7=0.  If  the  radar  axis 
is  stabilized  toward  S,  the  error  voltage,  averaged  over  a  long 
period,  also  vanishes.  The  aspect  angle  0  is  the  angle  between 
the  x~axis  and  the  LOS  from  the  tracker  T  to  S.  The  direction 
of  the  LOS  is  gjven  by  the  unit  vector  Uq,  the  direction  at  right 
angles  to  it  by  the  unit  vector  U^.  Note  that  TS  and  TSi  are 
almost  parallel  and  the  angle  is  very  small  under  "far-field 

conditions."  The  remaining  quantities  are  indicated  in  Figure  1, 
and  are  self-explanatory. 

-  7  - 
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(2)  Bistatic  Case 

This  case  refers  to  the  action  of  a  seni -active  honing 

interceotor.  The  target  configuration  is  the  same.  However,  the 

single  aspect  angle  ©  must  be  replaced  by  two  angles,  ©t  and 

©r.  The  transmitter  angle  ©t  is  the  angle  between  the  illuminator 

LOS  and  the  x-axis;  the  receiver  angle  0  is  the  angle  between 

r 

the  seeker  LOS  and  the  x-axis,  as  indicated  in  Figure  1.  In  the 
definition  of  the  error  angles  e^,  the  LOS  from  the  seeker  to 
the  mean  centroid  S  forms  the  basic  reference  line  in  the 
bistatic  case,  since  it  is  the  seeker  receiver  that  measures  the 
error  voltage. 


IV.  GENERAL  DEFINITIONS  AND  SYMBOLS 


In  the  analysis  of  radar  cross  sections  as  applied  to 
the  theory  of  target  detection,  the  following  quantities  play  an 
important  role: 

o 

a  is  the  radar  cross  section  of  a  body,  in  m  . 

The  RCS  is  4ir  times  the  re-radiated  electro¬ 
magnetic  power,  in  watts  per  steradian,  divided 
by  the  irradiated  power  per  unit  area,  in  watts 
per  v? . 

o^  is  the  instantaneous  RCS  of  a  body  computed  by 
the  method  of  "relative  phase."  Where  there  is 
no  confusion  likely,  we  shall  simply  omit  the 
subscript  "p". 

Oj  is  the  RCS  of  the  individual  component  scatterer. 

a  is  the  mean  RCS  (Ordinarily,  this  is  a  time 

average.  o’  *  o(t)  .  However,  o  might  also 
be  considered  as  a  parameter  process  with  the 
aspect  angle  Q  or  the  frequency  f  representing 
the  independent  parameter.  The  mean  then  is  taken 
with  respect  to  the  parameter.  ) 

to  *  a  -  a  is  the  instantaneous  deviation  of  the  RCS 

from  its  m>;an. 

— 5-1/2 

a  .  -  (to  )  is  the  standard  deviation  of  the  RCS. 
(Note:  Reference  [5]  uses  the  symbol  "S".) 


*In  the  bistatic  case  one  lias  o  »  a  {0^ ,  0r3  and  the  average  can 

be  defined  as  a  «  a  (o. (q),  0  (q) }  where  the  averaging 

t  r  av 

operation  is  extended  over  the  parameter  q. 


3.0 


is  the  maximum  possible  BCS  (taken  over  all 
possible  RF-phases  for  fixed  parameters) . 
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In  the  analysis  of  radar  target  glint  as  applied  to 
the  theory  of  target  tracking,  the  following  definitions  and 
symbols  are  important: 


(The  various  stations  of  the  target  are  indicated  by 
vectors  or  matrix  columns.  Vectors  are  denoted  by  capitals, 
their  components  by  lower-case  letters.  Conventional  matrix 
symbolism  is  used;  in  particular,  primes  designate  transposed 
matrices.) 


position  of  station 


X  "  [*»  y]  «  position  of  instantaneous  centroid  S 
(slow  AGC) . 

X  -  [x ,  y ]  -  position  of  mean  centroid  5  (slow  AGC) . 


AXj  -  Xa  -  X  -  deviation  of  station  SA  from  3 

AX  -  X  -  X  -  deviation  of  instantaneous  centroid 
from  mean. 

f 

(AX* AX  )  -  covariance  matrix  of  centroid  position 


% 


[cos  0,  -sin  0]  unit  vector  in  the  direction 
of  the  radar  LOS.  (0  -  line) 


[sin  0.  cos  0]  unit  vector  normal  to  the  9  -  line. 


r  -  range  from  radar  to  target,  in  units  consistent 
with  units  of  x. 


In  this  analysis  the  vectors  and  matrices  are  of  dimension 


"2” 


11 
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As  mentioned  In  Section  III,  It  is  convenient  to  refer  the  error 
angles  e  to  the  position  of  the  mean  centroid.  Hence 


(AX  -ty  n. 

e.  -  - - «  —  -  error  angle  of  station  S. 

1  r  r  l 


(see  Figure  2) 


(AX  -Gy) 

e  -  - - -  -  error  angle  of  instantaneous  centroid  (lb) 


and,  by  definition, 


e  -  0 


1/2 

,T 

esd  -  <e  >  • 


Simple  summations  extend  over  all  scatterers  Si  with 
1,  2,  ...,  N.  We  shall  use  the  following  abbreviations: 


£-? 


Double  summations  extend  over  all  combinations  of  i,  j  =*  1,  2, 
. . .  ,  N  except  terms  with  the  same  index: 


E  E  ~  E  -  2  l  l  =  2  i 

1  <  i  /  j  <  N  i/j  1  <  i  <  j  <  N  i  <  j 
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V.  SPECIAL  DEFINITIONS  AND  SYMBOLS  APPLICABLE 
TO  RICE'S  RANDOM  NOISE  THEORY 

For  the  application  of  S.O.  Rice's  "Random  Theory  of 
Noise"  (Reference  [10]),  a  few  statistical  definitions  are 
necessary. 

Let  z  -  z(q)  be  a  parameter  process.  (In  our  case, 
z  will  be  the  RCS  process  a( q) ,  while  the  independent  parameter 
q  could  be  the  time  t,  the  aspect  angle  O,  or  the  radio 
frequency  f.)  z  is  the  mean  of  z(q)  averaged  over  q,  and 
<z>  is  the  ensemble  mean  of  z(q).  For  homogeneous,  ergodic 
processes,  one  has: 

z  =  <z>.  (2) 

The  RCS  process  is  a  locally  homogeneous  process 
(locally  stationary  if  q  *  t) ,  and  equation  (2)  is  valid  at 
least  for  a  short  range  of  the  parameter  q.  Let 

Az  «=  z(q)  -  z  . 

Then  Az  is  an  unbiased  random  process. 

The  quantity  N  =  N{Az|Aq,  q,  p,...}  is  defined  as 
the  number  of  zero  crossings  with  positive  slope  of  a  sample 
process  Az(q)  in  an  interval  I  ranging  from  q  to  q  +  Aq. 

N  itself  is  a  random  variable  that  depends  on  Aq,  q,  and  possibly 
on  some  other  parameters,  such  as  p,  ...;  N  -  ITC  Az  |  Aq ;  q,  p,  ...} 
is  the  mean  of  N.  AN  -=  N  -  N  is  the  deviation  of  zero  counts 
of  a  particular  sample  process  from  its  mean  value. 

-  if,  Uz;  p,  ...}  -  p.  •••> 


* 


See  Appendix  A. II. 
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is  the  aean  count  of  zero  crossings  per  unit  of  q.  17  is  a  density 

-1  ,  ,  o  1/2 

with  diaension  q  .  lAqJ  -  (AN  )  is  the  standard 

deviation  of  counts  of  zero  crossings  with  positive  slope,  taken 
over  the  interval  I. 


A  properly  defined  standard  deviation  of  zero  crossing 
density  (or  rate  of  counts)  would  be: 


) 


q 


1/2 

f  lia  (ANUqjrl 
[Aq-0  Aq  ] 
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VI.  BASIC  ASSUMPTIONS  AND  PRINCIPLES  USED  IN  THE  ANALYSIS 
The  following  basic  assumptions  were  utilized  in  our 

analysis: 

Basic  Assumptions 

(1)  The  target  is  a  slender,  axially  symmetrical 

vehicle , 

(2)  it  consists  of  many  independent  scatterers,  that 
are  randomly  distributed  alongside  the  body, 

(3)  the  linear  dimensions  of  the  target  and  of  the 
distances  between  individual  scatterers  are  large  compared  to  the 
radio  wavelength  X, 

(4)  shadowing  effects  and  multiple  scattering  can  be 

ignored, 

(5)  the  range  r  from  the  radar  to  the  target  is 
large  compared  to  the  target  length, 

(6)  the  target  lies  within  the  regime  of  linear  error 
patterns  of  the  receiving  antenna, 

(7)  the  antenna  axis  is  space  stabilized  in  the 
direction  of  the  mean  target  centroid, 

(8)  the  receiver  uses  slow  AGC  and  square-law 

,  ** 
detection , 

(9)  receiver  noise  and  other  types  of  interference 
are  ignored. 

Assumptions  (1)  through  (4)  are  typical  for  the  special  kind  of 
targets  which  we  investigated,  assumptions  (5)  to  (9)  pertain  to 
the  properties  of  antenna  and  receiver  and  to  their  geometrical 
relationship  to  the  target.  All  these  premises  are  very  useful 
in  simplifying  the  analytical,  statistical  and  computational 
part  of  this  study,  especially  if  combined  with  a  few  mathematical 
and  physical  principles  which  will  be  listed  below  and  more  fully 
discussed  in  the  appendix. 

*  ______________________ 

See  discussion  at  the  end  of  Section  A. IV. 

cf.  remarks  in  Section  A.V.3. 
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Assumption  (1)  enables  us  to  reduce  the  problem  from 
a  3-dimensional  to  a  2-dimensional  job.  Assumption  (9)  implies 
that  receiver  noise  and  target  scintillations  are  statistically 
independent  of  one  another  and  therefore  can  be  treated  separately 
and  then  combined. 


Mathematical  and  Physical  Principles  which  Simplify  the  Analysis 


(a)  Far -Field  Scattering, 

(b)  Born -Approximation, 

(c)  Method  of  Random,  Relative  Phase, 

(d)  Central  Limit  Theorem, 

(e)  Principle  of  Local  Stationarity , 

(f)  Lord  Kelvin's  Principle  of  Stationary  Phase. 


Items  (a) ,  (to)  and  (c)  were  effectively  used  by  the 
Conductron  Corporation  (cf.  Ref.  [5])  and  will  be  discussed  in 
Appendix  A. I.  The  central  limit  theorem  (item  (u))  is  a  conse¬ 
quence  of  assumption  (2).  It  enables  us  to  utilize  Gaussian 
statistics.  This  property  and  the  principle  of  local  stationarity 
(item  (e))  entitle  us  to  make  use  of  some  results  of  Rice's  random 
noise  analysis.  Principle  (e)  was  introduced  by  us  for  this 
purpose.  It  is  a  consequence  of  assumptions  (2)  and  (?)  and  will 
be  described  in  Appendix  A. II.  Lord  Kelvin's  principle  of  stationary 
phase  (item  (f))  is  very  useful  in  simplifying  the  complexity  of  the 
numerical  analysis.  In  our  case,  some  of  the  final  expressions 
consist  of  fourfold  summations,  running  over  all  the  coordinates 
of  the  individual  scatterers.  If  we  deal  with  50-100  scatterers, 
the  number  of  terms  required  for  one  single  quantity  as  a  function 
of  a  set  of  fixed  parameters  (say  frequency,  aspect  angle  and 
polarization)  may  tax  the  capability  of  a  modern  electronic 
computer.  The  principle  of  stationary  phase  enables  us  to  reduce 
quadruple  sums  to  double  sums.  Details  of  this  method  will  be 
found  in  Appendix  A  TII. 
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VII.  SUMMARY  OF  MATHEMATICAL  RESULTS 

The  results  will  be  divided  into  two  main  sections  7.1 
and  7.2  entitled  "Simple  Statistics"  and  "Refined  Statistics" 
respectively.  The  former  contains  formulae  which  were  obtained 
previously  by  other  researchers  (cf.  Refs.  [5]  and  [19]),  the 
second  main  section  applies  the  methods  of  Rice's  analysis  to 
the  scintillation  problem.  The  results  are  believed  to  be  new, 
to  the  best  of  our  knowledge. 

Each  main  section  is  divided  into  two  subsections, 
e.g. ,  7.1.1  and  7.1.2.  The  former  deals  with  quantities  of  the 
detection  phase,  primarily  averages  of  radar  cross  sections,  the 
second  subsection  derives  parameters  which  are  useful  in  the 
tracking  problem.  We  are  concerned  with  the  erratic  meander 
of  the  target  centroid  in  the  radar's  field  of  view  (say  in 
azimuth  and  elevation) . 

Logically  a  third  and  fourth  subsection  should  follow 
the  second  one  and  produce  the  corresponding  parameters  in  range 
and  range  rate  (doppler) .  But  these  problems  were  not  part  of  our 
present  assignment.  They  will  be  analyzed  ±n  future  projects. 

(cf.  remarks  of  Section  VIII). 

7. 1  Simple  Statistics 

7.1.1  Radar  Cross  Sections  of  Fading  Targets 

Mean  Cross  Section 

a( 9)  =  £  a.  (0)  (3a) 

i  1 

In  the  bistatic  case,  the  quantity  0  .  is  replaced 
two  angles  0.  and  0  ,  as  in  equation  (3b)  and  all  following 

v  l 

formulae. 

V  -  £  ct.  (0t ,  0r)  (3b) 
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In  particular, 


Standard  Deviation  of  Apparent  Centroid  in  x -Direct ion 

i  r  r  —2} 

xsd  -  Vj<xt  -  XJ  - 2x)  J 


Standard  Deviation  of  Apparent  Centroid  in  y-Direction 
i  T  it1/'2 

y°<  ■  To  l  h  ,y‘ +  yj  - 2?)  J 


(8b) 


Standard  Deviation  of  Measured  Error  Angle  (Slow  AGO) 


» 

— <  J]  cr.cr.  (x.  +  x.  -  2x)  sin  © 

2vo  li/j  iJl  1  3  r 


+  (yt  +  yj  -  2y)  cos  ©r 


l  Refined  Statistics  (Density  of  Mulls) 

7.2.1  Radar  Cross  Sections  of  Fading  Targets 


Mean  Number  of  Lobes  per  Unit  of  Parameter 


N  Ua;q,p}  -  — 


[&  ^ 


(  f  y 


In  application  to  specific  parameters  and  cases, 
Table  1  serves  for  the  computation  of  Gq,  Table  2  for  the 
computation  of  the  term  qF^j . 
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VALDES  OF  G 

_ a 


Parameter  q 

Units 

G 

q 

Aspect  angle  0 

rad 

k  1 

2v  *  X 

deg 

k 

360 

Radio  frequency  f 

Hz 

1 

c 

MHz 

106 

c 

Time  t 

sec 

k 

b  8 

where 

k  -  ^  is  the  wave  number, 

A.  c  c 

c  -  speed  of  light , 

1 

0  -  rate  of  change  of  aspect  angle  in  radians/sec. 

The  units  of  X  and  c  must  be  compatible  with  those  of  a 
a  and  o. 
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Table  2 
VALUES  OF  F. 


Parameter  q 

1 

Monostatic  Case 

Bi static  Case 

0 

or 

t 

2[ (x±-Xj ) sin  0  +(yi-yj)cos  0] 

[(xi~Xj)(sin  ©t+sin  ©r)  + 

+  (yi-yJ)  (cos  ©t+cos  ©r>] 

f 

2[(xi-Xj)cos  0  -(y^yjjsin  0] 

[  (x .  -x.)(cos  O.+cos  0  )  - 
i  j  t  r 

-(y^yjMsin  ©t+sin  ©r)] 
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In  the  y -"Direction 


H,Uq;q.lj  -  2^-  [  £  °i°j  (y1+yj"2y)2(qFij)2] 


1/2 


(lib) 


While  the  last  two  quantities  are  not  directly  observable, 
the  following  quantity  can  be  measured  from  the  receiver: 

Mean  Number  of  Error  Fluctuations  per  Unit  of  Parameter 


!Lle;q,p} 


2ore 


sd 


■2x)  sin  0  + 
r 


U  °^h+xJ 

+  <y1+yj-2y)  cos  s,]2^.,)2  J 


1/2 


(12) 


One  example  follows: 


Mean  Number  of  Error  Fluctuations  per  MHz  at 
rixecfi  ligle  in  the  Monostatic  Case 


Wf{e;0} 


10' 


core 


sd 


+Xj-2jD  sin  0  + 


u 

+  (yA+yj-2y)  cos  o]2  • 

% 

•  C(xA-Xj)  cos  0  -  (y^-y^) sin  O]2  | 


(12a) 
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Application  of  Analysis 

The  Conductron  Corporation  of  Ann  Arbor,  Michigan, 

has  applied  the  mathematical  results  of  Section  VII  to  a  series 

$ 

of  specified  missiles.  A  set  of  numerical  calculations  was 
performed  for  the  output  quantities  listed  below  as  functions 
of  various  RF-wave lengths,  aspect  angles  and  polarizations: 


0 

eq. 

(3a) 

aMAX 

eq. 

(4) 

asd 

eq. 

(5) 

X 

eq. 

(6a) 

Xsd 

eq. 

(8a) 

NqUa;q,p] 

eq. 

(10) 

(1.) 

for 

q  -  0, 

p  -  0 

(2.) 

for 

q  -  t. 

p  -  0 

N  Ux;q,p} 

eq. 

(11a) 

(1.) 

for 

q  -  0, 

p  "  0 

(2.) 

for 

q  “  f. 

p  =*  0 

At  the  same  time,  the  Applied  Physics  Laboratory  under 
the  sponsorship  of  ARPA  carried  out  a  parallel  RCS  measurement 
program  at  Holloman  Air  Force  Base  by  static  ground  tests  on  models 
of  the  same  specified  missiles.  Again  various  RF -wavelengths , 
aspect  angles  and  polarizations  were  utilized.  What  conclusions 
can  be  reached  from  a  comparison  of  the  theoretical  and  experimental 
results? 


* 


Reference  [19],  (c)  and  (d) . 
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Before  answering  the  last  question,  a  few  general 
remarks  are  in  order.  Ve  deal  with  very  complex  and  extremely 
phase-sensitive  structures  in  the  microwave  regime.  At  first 
sight,  some  of  the  output  quantities  appear  to  fluctuate  wildly 
and  randomly  as  functions  of  certain  parameters  or  i  .  response 
to  slight  changes  of  the  configuration.  It  is  difficult  to  assign 
a  well-shaped  spectrum  to  the  output  signals.  A  flat  spectrum 
seems  to  be  the  closest  choice  for  them  at  best.  This  characteristic 
feature  is  confirmed  by  the  analysis  (see  Appendix  A. III).  It  is 
even  more  difficult  to  read,  extrapolate  or  guess  an  "instantaneous 
lobing  frequency"  from  the  inspection  of  the  test  records.  It 
was  indeed  this  difficulty  which  led  to  the  present  investigation. 

In  tbe  light  of  these  reaarks,  the  following  conclusions  can  be 
stated: 

(a)  There  exists  a  general  qualitative  agreement 
between  tbfe  computed  and  measured  quantities. 

0 b)  In  particular  theoretical  and  experimental 
results  show  the  same  trends.  For  instance, 
the  lobing  frequency  is  lowest  "head-on"  i.e., 
in  the  direction  where  the  aspect  angle  0  is 
zero,  it  rises  with  increasing  aspect  angle  and 
reaches  one  or  several  maxima  in  oblique  directions. 
The  maxima  are  not  fixed  for  the  same  missile  if 
the  radio  frequency  or  polarization  are  varied. 

As  expected  the  lobing  frequency  Ng  increases 
with  higher  radar  frequencies  f,  while  the  other 
parameters  are  held  fixed. 


The  discussion  of  the  class  of  deterministic  ana  random  functions 
which  alljw  an  "instantaneous  frequency"  will  be  taken  up  in  a 
future  report. 
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(c)  Provided  a  diligent  mathematical  and  numerical 

analysis  is  performed,  it  appears  that  the  accuracy 
of  the  important  output  quantities  is  sufficient 
for  a  statistical  radar-design  and  performance 
study  concerned  with  detection,  track  and 
identification  of  scintillating  targets. 


Recommendations 

The  computed  quantities  referring  to  the  meander  of 
the  centroid,  such  as  x  (eq.  6a),  xg<J  (eq.  8a),  N^{Ax;q,p] 

(eq.  11a)  et  al. ,  have  not  been  verified  as  yet  by  experimental 
measurements.  This  should  be  done,  preferably  in  dynamic  flight 
tests. 

As  mentioned  above  this  analysis  should  be  extended 
to  include  the  scintillation  of  the  outputs  of  the  range  and 
range  rate  (doppler)  circuits. 

The  investigation  should  be  generalized  to  incorporate 
targets  containing  non-random  features,  such  as  periodic  sections, 
uniformly  distributed  scatterer  segments,  etc.  (cf.  remarks  ^t 
the  end  of  Sec.  A. III). 

The  effects  of  body  vibrations  and. of  the  slower 
oscillations  of  the  stability,  control  and  guidance  loops  on 
the  radar  scintillation  should  be  examined  (cf.  discussion  in 
Sec.  A.V.l). 

Finally,  the  usefulness  of  special  devices  (RF  frequency 
modulators,  adaptive  narrow-band  filters  with  variable  center 
frequencies,  etc.)  should  be  studied.  The  objective  would  be  to 
recognize  the  size,  overall  shape  and  attitude  of  unknown  objects 
or  to  identify  the  signature  of  well  known  targets.  Some  of  these 
problems  are  being  studied  at  the  present  time  and  will  be 
documented  in  the  near  future. 
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A. I.  PRINCIPLES  USED  IN  THE  DERIVATION  OF 
THE  ANALYTICAL  RESULTS  OF  SECTION  7.1.1 

The  mathematical  results  of  Section  7.1  entitled 
"Simple  Statistics,"  Subsection  7.1.1  "Radar  Cross  Sections 
of  Fading  Targets"  were  derived  by  the  research  team  of  the 
Conductron  Corporat/on.  For  details  of  this  analysis,  Refs. [5] 
and  [19]  should  be  consulted.  Here  we  shall  discuss  the  basic 
principles  utilised  in  this  work. 

(a)  Far-Field  Scattering, 

(b)  Born  Approximation, 

(c)  Method  of  Relative,  Random  Phase. 

Principle  (a)  is  based  on  our  assumptions  (5)  and  (3) 
of  Section  VI.  Tt  simplifies  the  problems  of  wave  physics  and 
in  many  cases  allows  the  use  of  plane  wave  fronts. 

Principle  (b)  is  well  known  from  scattering  problems  in 
quantum  mechanics.  There  it  becomes  applicable,  if  the  interaction 
coefficients  are  so  small  that  double  ana  multiple  scattering 
involving  squares  and  higher  powers  of  the  interaction  constants 
can  be  ignored.  In  our  case  multiple  scattering  can  be  omitted 
since  we  deal  with  bodies  of  simple  shapes  (assumptions  (1)  and 
(4)  of  Section  VI). 

Principle  (c)  or  the  method  of  relative  and  random 
phase  will  be  described  next.  Following  the  analysis  of  Ref.  [5], 
page  25  a.f.  and  using  the  notation  of  Section  V,  one  obtains  for 
the  instantaneous  RCS  of  the  target: 

N 

Op  -  I  7  toj)172  exp<l*j>|2  (A.I.l) 
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where  in  the  wonostatic  case 

-  -2k[(x^-x)cos  0  -(y^-y)sin  0]  +  $ 
in  the  bistatic  case 


-  -k[  (Xj  -x)  (cos  ©t+cos  Or)-(y^-y)  (sin  Ot+sin  ©r>]+  $ 

(A. I, 2b) 

k  is  the  wave  number  explained  at  the  bottom  of 
Table  1. 

$  is  an  absolute,  instantaneous  phase  of  the  mean 
centroid. 

The  most  difficult  part  of  the  analysis  is  the  job  of  finding 
the  magnitude  of  the  individual  scattering  cross  sections  o^. 

The  derivation  cf  ap  and  a  proceeds  in  three  steps.  To 
describe  the  method,  we  follow  Ref.  T5]  closely  and  verbatim  in 
BRny  instances: 

First  Step  -  Modeling 

We  consider  the  body  under  investigation  as  an  assembly 
of  many  components  each  of  which  can  be  geometrically  approximated 
by  a  "simple  shape"  in  such  jl  way  that  the  RCS  of  the  simple  shape 
approximates  the  RCS  of  the  component  it  models. 

Second  Step  -  Computation  of  Component  RCS 

It  is  assumed  that  the  RCS  of  the  "simple  shape"  can 
be  computed.  The  rigorous  method  would  require  a  solution  of 
Maxwell '8  equations  with  known  but  complicated  boundary  conditions. 

If  an  exact  solution  is  not  known,  one  must  find  an  approximate 
solution  by  some  special  method. 
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Third  Step  -  Combination  of  Individual  Component  RCS's 


The  final  step  involves  the  proper  combination  of  the 
’’component  cross  sections"  to  yield  the  estimate  of  the  cress 
section  of  the  entire  body.  Here  again  it  is  often  necessary  to 
find  a  simplifying  mathematical  scheme  such  as  the  method  of 
"relative  random  phase."  It  must  be  emphasized  that  mathematical 
approximations  must  be  Justified  on  physical  grounds. 


Rewriting  eq.  (A.I.l)  one  gets  in  the  monostatic  case: 
N  N 

ap  *  /L^  Z  (a±opl/2cos  2k[  (Xj  -x^cos  Q  -(yj-yi)sin  ©] 


In  the  notation  of  Section  IV: 

a  -  ]T  a.+  2  T  (a.a,)1//2cos  2k[(x.-x.)cos  ©-(v.-y.)sin  0]  (A. 1.3) 

p  i  i<j  J  J  31 

It  is  seen  that  the  absolute,  instantaneous  phase  of  the 
centroid  disappears  from  the  equation  and  therefore  becomes 
irrelevant.  This  fact  justifies  the  term  "relative  phase." 

A  particular  procedure  of  combining  individual  component  cross 
sections  is  called  the  "method  of  random  phase."  It  is  based  on 
the  assumption  (see  p.  26  of  Ref.  [5])  that  the  many  different  ^ 
are  randomly  distributed;  then  upon  averaging  over  one  obtains 

the  expression 

a  ~  £  a,  (A. 1.4) 

i 

which  is  identical  with  eq.  (3a)  of  Section  VII. 

The  reasoning  which  led  to  eq.  (A. 1.4)  appears  some¬ 
what  vague.  However,  one  has  the  feeling  that  the  result  is 
physically  correct.  One  might  try  to  strengthen  the  mathematical 
argument  by 
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(a)  introducing  statistical  distributions  for  the  x^ 
(say  Poisson  distributions)  and  independent  but 
identical  distributions  tor  the  (say  Rayleigh 
distributions) , 

(b)  going  to  the  limit  -  oo  ,  (L  is  the  length  of 
the  target) 

(c)  replacing  sums  by  integrals  and  averaging. 


Then  one  vinds  up  with  expressions  containing  Dirac's 
delta  functions,  giving  rise  to  analytical  difficulties. 

An  alternate  method  of  combining  wildly  fluctuating 
signal  components  will  be  discussed  in  Section  A. II. 
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A. II.  PRINCIPLES  AND  METHOD®  USED  IN  THE 
ANALYSIS  OF  SECTION  7.2 


Referring  back  to  the  list  of  principles  enumerated 
in  Section  VI,  we  shall  make  use  of  the  following  items: 

(d)  Central  Limit  Theorem, 

(e)  Principle  of  Local  Stationarity , 

(f)  Lord  Kelvin's  Principle  of  Stationary  Phase  and 
its  Extension, 

Item  (d)  is  based  on  assumption  (2)  of  Section  VI.  If 
we  deal  with  many  independent  scatterers,  whose  individual  RCS's 
follow  independent  but  equal  distributions,  and  if  some  other 
conditions  are  met  (such  as  the  existence  of  first  and  second 
moments)  ,  the  RCS  of  the  combined  signal  in  the  limit  N  -*  ® 
will  approximate  a  Gaussian  distribution.  This  property  is 
important  because  it  enables  us  to  make  use  of  seme  results 
obtained  by  S.  0.  Rice  in  the  field  of  the  statistics  of  the 
count  of  zero  crossings. 

Item  (e)  is  a  new  concept  which  was  introduced  here 
especially  for  the  purpose  of  coping  with  the  type  of  signals 
encountered  in  radar  target  echoes.  These  signals,  in  the  long 
run,  or  averaged  over  a  long  interval  of  time,  display  a  broad, 
flat  spectrum.  .  Yet  we  know  these  signals  bear  an  instantaneous 
statistical  character  which  changes  slowly  with  the  parameter 
(aspect  angle  9,  time  t  or  frequency  f) . 

The  concept  of  "local  stationarity"  shall  be  explained 
on  hand  of  the  signal  ap  (cf.  eq.  A.I.l). 

N 

<7p  -  (aj^1/2  exP(1$j>|2  <A 


* 


Reference 


[14],  p.  215. 


.11.1) 
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For  simplicity  we  choose  the  monostatic  case  and  a  one- 
dimensional  distribution  of  scatterers.  Accordingly: 


^  -  -2k[ (Xj -x)  cos  0]  +  * 


(A. II. 2) 


Then  (A.II.l)  can  be  rewritten  as  follows: 


ap(0)  -  ^  (a^a^)^2  cos  [2k(Xj-xi)  cos  0] 


(A. II. 3) 


0  is  the  independent  parameter  of  the  real  signal 
process.  Now  we  introduce  two  additional  independent  virtual 
variables  t  and  6: 


a  (9,1,6)  -  £  (o.< J.)1//2  cos[2k(x.-x.)  cos (9  +  t  +  6)] 

P  i^j  1  3  J  1 


(A. II. 4) 


t  is  the  parameter-shift  variable  in  a  statistical  ensemble 
with  the  range  0  <  t  <  2*  and  the  probability  density: 


p(T)  -  £ 


(A.II.5) 


6  is  the  parameter  in  a  locally  stationary  ensemble  with  the 
range  -A  <6  <  + A  and  the  probability  density 


"  2A 


(A. II. 6) 


Irrespective  of  bow  small  we  choose  A,  we  can  always 
find  a  large  enough  k  or  k  *  (x^-x^),  so  that  the  argument  or 
phase  of  every  cos-term  describes  the  full  circle  from  0  to  2ir, 

even  several  times.  In  the  limit,  if  k(Xj-x^)  -•  oo  ,  the  process 
a  (6)  appears  as  a  full  random  process  covering  infinite  phase 

Mr 

shift  as  6  goes  through  its  range  -A  <  6  <  +A.  In  the  following, 
we  shall  always  assume  that  we  are  close  to  this  idealized  situation. 
For  lack  of  other  suitable  terms,  we  shall  call  a  process  of  this 
type  a  "peculiar  process,"  i.e.  ,  peculiar  to  certain  types  of 
radar  echoes. 
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Indeed,  the  peculiar  processes  share  with  the  flat, 
white  noise  processes  the  broad  spectrum.  However,  they  appear 
to  be  more  closely  related  to,  but  not  identical  with  the  noise 
processes  of  Beference  [16],  eq.  (12). 

The  peculiar  processes  are  not  stationary  over  the  full 
ranges  of  their  parameters  (0,  t) .  However,  for  fixed  0  and  t, 
we  may  define  a  "local  stationarity , "  by  taking  a  "microscopic 
look"  so  to  speak,  at  the  small  range  -A  <  5  <  +A  of  the  local 
ensemble.  By  making  a  transformation  of  the  form 


we  obtain  a  process  of  the  parameter  6  which  is  stationary  in 
the  range 


~v  <  8  <  +w. 

This  procedure  is  similar  to  Newton's  idea  of  forming  derivatives. 
To  any  point  Px  of  a  curvilinear  function  choose  an  infinite¬ 
simally  near  point  P2  on  the  curve.  Connect  Pj^  with  P2  by 
a  straight  line.  The  infinitesimal  segment  PjPg  can  then  be 
blown  up  by  extending  it  to  a  tangential  line.  It  is  well  known 
that  derivatives  do  not  exist  for  certain  classes  of  real,  even 
continuous  functions.  In  the  same  manner,  stochastic  processes 
do  not  necessarily  possess  local  stationarity,  but  the  peculiar 
processes  enjoy  this  property. 

The  virtual  variable  t  is  used  in  connection  with 
ordinary  stationary,  statistical  quantities,  such  as  the  auto¬ 
correlation  function,  spectrum,  etc.  The  virtual  variable  6 
is  used  in  connection  with  locally  stationary  properties,  such 
as  "instantaneous  frequency,"  "average  number  of  zero  crossings," 
etc. 

We  can  speak  of  several  types  of  statistical  averages 
for  a  quantity  z  of  such  a  process: 
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The  parameter  ttime)  average  is  indicated  by  z. 
The  ensemble  average  is  denoted  by  <  z  >.  It  is  veil  known 
that  for  stationary,  ergodic  processes  one  has: 


s  «  <  z  > 


(A. II. 7) 


Ve  have  aany  real  coaponent  scatterers.  Therefore  we 
could  also  speak  of  an  assembly  average  z.  In  otr  vork  assembly 
averages  are  not  required,  but  total  suas  over  the  assembly  are, 
which  in  the  simplest  cases  become  tfz.  For  instance  one  gets: 

✓  oj’-  H|<j^1/2  exp  <i*j)| 


The  local  ensemble  average  {z(0^(t>} 

+A 


{z(0,T>} 


lim  1 
A-0  2l 


dd  z(0  +  t  + 


fi) 


is  defined  by 


(A. II. 8) 


It  is  important  to  note  that  for  our  type  of  processes 
the  assembly  sums  change  with  finite  shifts  of  the^real  parameter 
0  and  of  the  ensemble  parameter  t.  A  rigorous  mathematical 
analysis  would  have  to  vork  around  Dirac's  d -functions,  since 
the  processes  considered  have  a  flat  spectrum  in  the  res’  and 
ensemble  domain.  Our  present  objective  is  limited  to  reducing 
expressions  which  require  exorbitant  calculational  work ) to  a 
tractable  computational  fora.  Therefore  we  shall  omit  further 
mathematical  and  theoretical  exercises  at  this  time. 


Generalization  of  the  Principle  of  Stationary  Phase 

The  principle  of  stationary  phase  was  originally 
enunciated  by  Lord  Kelvin  as  a  mathematical  artifice  to  approxi¬ 
mate  infinite  integrals  of  wildly  fluctuating  functions  of  the 
type  occurring  in  the  hydrodynamic  theory  of  impulsive  disturbances 
or  wave  groups  of  small  spatial  extension.  These  waves  possess  a 
broad  spectrum.  This  principle  was  later  cast  in  a  mathematical 
rigorous  form  (see  Ref.  [15]). 
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We  shall  state  this  principle  in  2  slightly  more 
general  forms  applicable  to  our  locally  stationary  processes. 

First  Version  of  Kelvin's  Principle 

Given  a  finite  function  f{*}  of  bounded  variation, 
such  as  the  cos-function.  The  argument  is  called  phase.  The 
latter  itself  is  a  function  of  two  sets  of  variables:  {m}  and 
{x,y , . . .  } 


*  -  •*  (m;x  ,y  , , . .) 


m  is  continuous  and  x,  y,  ...  assume  discrete  values  only. 

We  form  infinite  integrals  over  the  first  variable 
and  finite,  but  large  summations  over  the  second  set. 


(A.II.9) 


If  the  phase  $  fluctuates  wildly  as  a  function  of  m, 
then  one  can  obtain  a  good  approximation  of  .  E 

(a)  either  by  integrating  over  a  short  range  of  m: 

nr  -  A  <  m  <  m  +  A  and,  or 
o  —  —  o 

(b)  by  summing  over  a  smaller  subset  x,  ,  x0  , . . . ,  y,  ,  y„  , . . . 

o  o  Ao  ^o 

of  the  second  set,  jo  long  as  the  stationarity  of  phase  holds  for 
the  two  subsets  (a)  and  (b) ,  i.e,  ,  so  long  as  the  following 
equation  is  valid: 

“  {mQ;  x^  ,y  ,...}=  0  for  mQ  and  all  x±  ,  y.  .  (A. II. 10) 

oJo  o  Jo 

Then 
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E"J  da  £  ]T  f{<Km;  x^y^...)}  (A. II. 11) 

a  -A  i“io’ 
o 

The  second  part  (b)  of  the  proposition  constitutes  the 
generalization  of  Kelvin's  principle. 

Example :  Let  us  apply  principle  (f)  to  our  locally  stationary 
signal  Op  (eq.  A. II. 4).  The  dummy  variable  <5  is  substituted 
for  a  and  Xj,  for  xA,  y^.  The  variables  0  and  t  are 
not  suaaed  here.  In  our  case  6  runs  through  a  very  short  range 
already.  Therefore  part  (b)  is  expected  to  play  a  major  role  in 
the  siaplification  process. 

Equation  (A. II. 10)  becomes 


0  «  ^  [2k(xA  -x.  )  cos  (0  +  t  +  6q)]  or 
o  Jo 


(xA  -x.  )  sin  (0  +  t  +  6q) 
o  Jo 


(A. II. 12) 


Since  0  and  t  are  arbitrary  this  equation  can  only  be  satisfied 
if  iQ  -  jQ  or  starting  with  the  original  set  £ i , j  5  one  forms 
the  subset 


The  local  ensemble  average  of 


a  becomes 
P 


(A. II. 13) 


Since  the  parameter  0  disappeared  from  the  right  band  side  of 
eq.  (A. II. 13),  this  expression  is  also  the  parameter-average  o: 


(A. II. 14) 
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This  is  the  result  of  Ref.  [5]  and  the  content  of  eq. 

(A. I. 4).  In  later  examples  we  shall  form  local  ensemble  averages 
that  not  identical  with  the  parameter  averages.  The  use  of 

the  extension  of  Kelvin* s  principle  somewhat  strengthens  the 
intuitive  argument  given  on  pg.  26  of  Ref.  [5],  The  validity  of 
Kelvin's  principle  was  proved  for  a  large  class  of  oscillating 
functions  by  G.  N.  Watson.*  The  extension  of  this  principle  should 
be  treated  in  a  more  rigorous  form  too.  But  this  job  exceeds  the 
scope  of  our  present  assignment. 


Second  Version  of  Kelvin's  Principle 


We  shall  cow  describe  a  second  version  of  Kelvin's 
principle  of  stationary  phase  that  will  be  utilized  in  the  following 
section.  In  the  previous  case  the  phase  function  contained  2  sets 
of  variables:  a  continuous  variable  and  a  set  of  discrete  variables. 
Accordingly  the  total  summation  consisted  of  an  integration  and  a 
finite  summation  over  discrete  variables.  In  the  sequel  the  phase 
function  contains  two  sets  of  discrete  variables  and  the  total 
summation  consists  of  two  finite  summations  over  discrete  variables. 
In  particular  each  summation  could  be  a  double  sum  over  discrete 
variables,  as  follows: 


E  = 


f 


(A. II. 15) 


Again,  if  the  phase  9  fluctuates  wildly  as  a  function 

of  the  variables,  one  may  get  a  fair  approximation  of  E,  if  one 

sums  only  over  a  subset  {x.  ,x,  :  x  ,x  }  of  the  set  {x. ,x. ;x  ,x  ] 

io  j  o  no  mo  ^  i  j  n  m 

for  which  the  phase  is  stationary  or  constant: 

®(x.  ,x.  ;  x  ,x  )  =  constant  for  any  choice  of  (A. II. 16) 
io  jo  no’  mo 

the  x-variables  within  the  subset. 

It  turns  out  that  in  the  special  applications  of  the 
following  section,  the  constant  vanishes  due  to  the  fact  that  the 
component  scatterers  x^  are  assumed  to  be  distributed  in  a  random 
fashion,  i.e.,  umerically  placed  at  incommensurate  distances. 


* 


Reference 


[17] 
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A.  III.  APPLICATION  OF  RICE'S  THEORY  TO  THE 
COMPUTATION  OF  RADAR-CROSS  SECTION  NULLS 


The  radar  cross  section  a of  a  target  which  satisfies 
the  conditions  of  Section  VI,  can  be  computed  by  the  method  of 
relative  phase  and  is  given  by  eq.  (A. I. 3).  With  slightly  modified 
notation  it  becomes: 


Jp  -  +  £  (°icj)1//2  cos  2k  [Otj  -x^cos  0  -  (yj-y^sin  ©] 


(A. 1.3) 


Op  -  cip(©)  can  *>e  considered  as  a  stochastic  parameter  process 
with  the  aspect  angle  0  playing  the  role  of  the  independent 
parameter.  According  to  eq.  (A.  1. 4),  the  mean  cf  of  op  is 
given  by 


Therefore 


AO  -  a  -  a 

P 


1//2  cos  2k  [  (Xj  -x^)  cos  O 


(A. 1.4) 


(yj -yi)sin  O] 


(A. III. 1) 


is  an  unbiased  parameter  process.  Moreover,  it  approximates  an 
unbiased  Gaussian  process  for  the  following  reasons:  If  we  hold 
i  fixed,  the  r.h.s.  of  eq.  (A.III.l)  is  a  sum  of  N  terms  T^ 
for  j  -  1,2,  ...,  N.  Each  of  these  terms  follows  an  independent 
distribution  (according  to  assumption  2  of  Section  VI) .  Further¬ 
more  all  finite  moments  of  each  term  exist,  since  we  deal  with  a 
finite  number  of  scatterers  and  the  total  extension  of  the  target 
is  bounded.  Therefore  the  conditions  of  the  central  limit  theorem 
are  met  (see  Ref.  [14],  pg.  215).  Hence  each  term  TA  with 


approximates  a  Gaussian  distribution  as  N  grows  toward  infinity. 


Again  the  monostatic  case  is  analyzed  here  for  simplicity. 
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The  parameter  process  ha  can  be  written  as 

40  -  L  t* 

i.e.  ,  as  a  sum  of  N  Gaussian  processes.  Therefore  ha  itself  — 
according  to  a  well  known  theorem  of  statistics  —  will  be  Gaussian 
or  —  more  precisely  —  vill  approximate  a  Gaussian  process,  as  N 
tends  toward  infinity,  ha  is  not  a  stationary  process,  that  is, 
its  statistical  properties  are  not  invariant  toward  a  finite  shift 
of  the  parameter  0.  However,  we  learned  in  Section  A. II  that  it 
could  be  considered  as  a  locally  stationary  process.  For  any  fixed 
choice  of  0  we  can  form  a  local  ensemble  with  the  parameter  6 
and  range  -A  <  6  <  >A  for  which  ha  behaves  like  a  stationary 
process. 

The  process  Aa(0)  defined  by  eq.  (A.III.l)  and  described 
in  the  previous  disc»_c  *ion  shall  be  called  an  "unbiased  peculiar 
process."  To  bring  out  the  essential  features  of  the  analysis  and 
save  space  at  the  same  tire,  we  shall  assume  that  v ’*  r.h?  v  'r.poaent 
scatterers  are  distributed  over  a  linear  segment  The  equations  of 
Section  VII  (Summary  of  Mathematical  Results) ,  however ,  cover  the 
general  two-dimensional  case. 

The  Auto -Correlation  Function  ui  an  Unbiased 
Peculiar  Process 

The  second  version  of  Lord  Kelvin's  principle  of  stationary 
phase  (see  Section  A. II)  shall  now  be  applied  to  the  computation  of 
the  auto-correlation  function  p( t,0)  of  the  process  Aa(O) , 
c.f.  eq.  (A.III.l).  Using  the  average  over  a  local  ensemble, 
one  gets: 

+ A 

PAa( t,0)  -  2£  f  d<5  Aff(Otfi)  Aa(0  +  x  +  6)  (A. III. 2) 

-A 

Substituting  eq.  (A.III.l)  in  (A. III. 2)  and  reducing 
the  problem  to  one  dimension,  one  has: 
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or 


1/2 


x  cos [2k  cob  (ftffi) (x.hs.) ]cosf2k  cos  (0  +  6  +  t) (x  -x  ) ] 

j  1  an 

*10<T-8>  -  srj  <“  *1  nI  <v//2  <v//! 

-A  J 

x  (cos[2k  cos  (046)  (x,~x/)  +  2k  cos  (6  +  6  +  t) (x  -x  ) ] 

j  1  m  n 

4  cos[2k  cos  (046) (x.-x.)  -  2k  cos  (©  +  6  +  t)  (x  -x  )]} 

j  i  m  n 


(A. III. 3) 


Hie  radar  cross  sections  of  the  component  scatterers  satisfy 

&  certain  distribution.  It  is  well  known  that  the  strong  scintil¬ 
lation  of  radar  echoes  is  caused  primarily  by  the  interference  of 

phases  between  the  component  scatterers  rather  than  by  the  fluctuation 

♦ 

of  the  o^  of  the  individual  scatterers.  Be  it  as  it  may,  let  us 
assume  temporarily  that  the  are  constant  and  equal  and  their 

variations  are  absorbed  by  slight  changes  of  the  large  kxA 
quantities. 

Since  j  >  i  and  m  >  n,  the  coefficients  of  the  cos 
terms  in  the  square  bracket  of  the  second  line  of  eq.  (A. III. 3) 
are  positive  and  therefore  the  whole  expression  in  the  second 
line  will  vanish  since  the  heavily  fluctuating  terms  cancel  one 
another.  A*  similar  situation  prevails  for  the  expression  of  the 


This  property  is  discussed  in  Hef.  [1],  pg.  202.  The  fact  that 
the  individual  do  not  depend  on  0  is  explained  by  the  large 

width  of  the  scattering  lobe  patterns  of  the  individual  scatterers, 
in  contrast  to  the  narrow  lobes  produced  by  an  array  of  length  L. 
The  constancy  of  the  individual  04  as  a  function  of  the  radio 
frequency  is  pointed  out  in  Ref.  [22],  pg.  573. 
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third  line  so  long  as  t  >  0,  since  we  assumed  a  completely  random 
and  incommensurate  distribution  of  the  x^.  However,  if  t  »  0, 
the  situation  is  different.  Then  there  exists  a  subset 

txio'  xJo:  xno'  x.o!  of  the  8et  txf  xj:  V  X«J  for  whlch  tbe 
phase  (term  in  square  bracket)  is  constant,  even  zero  for  all 

combinations.  The  subset  is: 


i,  j  ■  1(  2 . N 


(A. III. 4) 


The  auto-correlation  function  becomes: 


P40(T'8)  * 


T  o4o.  for  t  -  0 
1  J 


(A. III. 5) 


for  t  /  0 


As  N  tends  toward  infinity,  p  (t, 0)  approximates 

ijO 

a  Dirac  delta -function  for  all  values  of  0.  Therefore  the 

spectrum  becomes  wide-banded  and  flat.  This  property  can  also  be 

£ 

proved  by  other  methods  and  certainly  was  expected.  However,  it 
is  highly  surprising  that  one  can  define  an  instantaneous  frequency 
or  rate  of  lobing  for  the  "peculiar  process."  This  concept  is  based 
on  the  noise  theory  of  S.  0.  Rice  (Ref.  [10]). 

The  Average  Rate  of  Nulls 

For  a  Gaussian,  stationary  unbiased  process  z(t)  with 
correlation  function  p  (t) ,  Rice  derived  the  quantity  N.{z),  i.e., 
the  average  number  of  zero  crossings  with  positive  slope  (number 
of  lobes)  per  unit  time: 


Compare  discussion  in  Appendix  A.V.2. 
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srt<*> 


1  , 
57  < 


PZT7T 


1/2 
^  T-0 


This  concept  can  be  carried  over  directly  to  Gaussian, 
locally  stationary,  parameter  processes.  For  the  unbiased  peculiar 
processes  one  obtains: 


(A. III.  6a) 


Mote,:  In  the  non-st atlonary ,  though  "locally  stationary"  processes 
all  quantities  depend  on  the  instantaneous  value  of  the  parameter 
».  The  surprising  property  of  the  peculiar  processes  mentioned 
above  can  now  be  explained  by  the  fact  that  while  p(t-*0)  and 
P(t-«0)  tend  toward  infinity,  their  ratio  remains  finite  and  gives 
the  desired  value. 

We  start  with  eq.  (A.  III.  3)  and  simplify  tlie  expression 
on  the  r.h.s.  somewhat  by  omitting  the  second  line  which  oscillates 
heavily  and  averages  to  zero  and  skipping  the  integral  over  the 
range  of  the  local  ensemble  variable  6  which  for  the  specific 
conditions  of  our  example  cancels  out. 

*,40<t’B)  -  2  X  E  (<,iaj)1/2  (Vm)1/2 


X  cos [2k  cos  0  (x.-x.)  -  2k  cos  (Of t)  (x  -x  ) ]  (A. III. 3a) 

j  i  in  n 

Differentiating  eq.  (A. III. 3a)  twice  with  respect  to  t,  gives: 

"  &  kW”  <v")1/2 

X  /cos  2k[cos  0(x  ,  -x. )  -cos (Of  t)  (x  ~x  ) ] • [2ksin (Of t)  (x  -x  )  3^ 
L  j  1  m  n  m  n 

f  sin  2k[cos  Q(x . -x . ) -cos (Of  t) (x  -x  ) ] • [2kcos(Of t)  (x  -x  ) ] > 

ji  mn  mnj 


(A. III. 7) 
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Nov,  we  keep  n  and  m  fixed  and  vary  i  and  j, 

Vith  our  special  conditions  (very  large  k,  random  distribution 
of  the  x^,  relative  constancy  of  oi)  the  summation  over  i,  J 
of  the  terms  in  the  second  line  which  fluctuate  wildly,  average 
to  zero  except  for  a  subset  (A)  {i0,  jQ;  nQ,  mQ)  of  £ i , J ;  n,m}: 


no 


mo 


io 


x 

■ 

X 


no 


mo 


(A)  n,  m-1,  2,  ...,N 


(A. III. 8) 


For  the  subset  the  phases  of  the  cos  terms  in  the  second 

line  disappear,  as  t  goes  to  zero.  The  same  reasoning  holds  for 

the  third  line.  Therefore  the  third  line  vanishes  since  sin  ♦  *  0 

n,m 

if  '♦  »  0.  Combining  these  results  and  interchanging  i,  j 

with  n,  m  one  finally  gets: 


NeUa;0) 


2v 


X^(/aJ/aJ)2[2k  sin  0(xj-xj)T 

$,(/w4 


1/2 


(A. III. 9) 


This  result  can  be  easily  extended  to  two  dimensions  by 

assuming  a  two-dimensional  distribution  of  the  scatterers,  random 

and  independent  in  x  and  y.  Furthermore  instead  of  using  the 

aspect  angle  0  as  the  independent  parameter,  one  may  substitute 

the  radio  frequency  f  or  the  time  t  or  any  ether  suitable 

variable  q.  In  the  auto-correlation  function  p.  ( t)  one  has 

to 

to  substitute  t  -  ,  t  -  Tt  or  t  -  for  t  -  tq.  Finally 

in  the  derivation  of  N  (Aa;q)  one  has  to  apply  the  second  derivative 
of  p( t  )  with  respect  to  t  .  The  results  are  summarized  in 

q  q 

eq.  (10)  combined  with  Tables  1  and  2  of  Section  7.2.  Note:  While 
numerator  and  denominator  of  the  r.h.s.  of  eq.  A. III. 9  tend  toward 
infinity  as  N  grows  over  all  bounds,  their  ratio  remains  finite 
and  yields  a  finite  and  reasonable  value,  as  was  shown  by  a 
numerical  analysis  carried  out  by  the  Conductron  Corporation. 
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Effect  of  Non-Random  Distribution  of  the  Scatterers 


If  the  scatterers  are  not  randomly  distributed,  ther 
the  reasoning  that  led  from  eq.  A. Ill, 3  to  eq.  A, III. 5, break*  loan. 
Supposing  the  scatterers  display  some  periodic  features  as  follows. 


ai+lp  “  °i 

L 

Xi+(A+l)p  ~  Xi+Ap  "  Xl+p  “  Xi  "  N+l 


where  p  is  a  fixed  integer  counting  the  number  of  scatterers 
within  a  period  and  A  is  a  running  integer  0  <  A  <  M  counting 
the  number  of  periodic  sections.  Then  one  can  enlarge  the  subset 
(A. III. 4)  by  adding: 


xno  "  xio  +  Ap 
*mo  ”  Xjo  +  ip 


(B) 


A  -  1,  2,  . . . ,  M 


Furthermore  there  will  exist  certain  discrete  values  of 
t,  say  Tg,  ....  i ^ ,  ...  for  which  additional  subsets 

Sp  S2,  ...,  S,,  ...  with  stationary  phase  are  associated,  for 
instance  for  one  can  construct  a  subset  3^,  as  follows: 


xio 

Xjo 


XJ 

xio 

Xjo  +  £p 


(S^)  i,  j  »  1,  2,  ...,  N 
l  is  fixed. 


(A. III. 10) 


* 

Each  index  k  appearing  in  the  subscripts,  should  satisfy  the 
condition:  0  <  k  <  N.  If  this  condition  is  not  met,  replace  k 
by  an  integer  k'  which  satisfies  this  inequality  and  the  following 
relation: 

k  *  —  k  (mod  N) 
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ao  that 

CC3  0  •  (x-x.)  -  cos  (e+T.)  (x  -x  )  -  0,  (A. III.  11) 

JO  lO  X  BO  DO 

In  addition  one  can  also  form  subsets  of  type  (B) . 

The  correlation  function  becones  a  sum  of  Dirac  delta 

functions: 


If 

P<T>  “  Iba<T'Ti)'‘1 


(A. III. 12) 


and  the  radar  echo  pattern  assumes  the  well  known  finger- like 
or  fan-like  appearance. 
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A.  IV.  EFFECT  OF  AGC  ON  THE  CHARACTER  OF 
CENTROID  OSCILLATIONS 

1.  Instantaneous  Centroid  for  Slow  AGC 

Of  — 

Let  k  -  -s—  be  the  wave  number  and  let  w.  -  k(P.~P  ) 

A  X  X  O 

for  i  -  1,  2,  N  be  the  relative  phase  shift  between  scatterers 

i  and  0.  PQ,  Pj,  P^,  Pjj  denote  the  two-^ray  paths  from 

tb*  radar  to  the  individual  scatterers  SQ,  S . ..,  Sj,  SN 

and  back  to  the  radar.  In  the  monostatic  case,  the  paths  run  from 
the  tracking  radar  to  the  scatterer  and  back  to  the  same  radar.  In 
the  bistatic  case,  the  paths  extend,  e.g. ,  from  the  illuminator  to 
the  scatterer  and  back  to  the  seeker  of  a  seaii -actively  homing 
missile.  SQ  is  a  reference  station,  say  the  nose  tip  of  the 
target  or  the  average  centroid  position. 

In  order  to  determine  the  position  of  the  centroid,  one 
has  to  define  two  vectors  V  and  V  and  one  square  matrix  M. 

A 

The  matrix  M  is  defined  by  its  general  element: 

-  cos  i,  J  -  1,  2,  ...,  N  (A.IV.l) 

where  4^  is  the  phase  angle  defined  by  eq.  (A. I. 2a). 

The  vectors  V  and  Vx  are  defined  by  the  following 

columns: 


/  0^ 

/ajx1 

^2*2 

V  - 

• 

• 

/3T 

• 

V  - 

X 

• 

• 

/ojx 

• 

• 

• 

• 

* 
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a2,  ...,  aA . aN  are  the  radar  cross  sections 

of  the  individual  component  scatterers  as  defined  in  Section  IV. 

The  quantities  x  and  are  used  in  the  present  section  as 

generic  terms  that  might  represent  either  the  x-  positions,  the 
y-  positions  or  the  error  angles  e  of  individual  scatterers 
and  centroids,  as  defined  in  Section  IV. 

Making  use  of  the  assumptions  listed  in  the  beginning 

of  Section  VI,  one  may  write  the  position  of  the  instantaneous 

centroid  x  for  slow  AGC  as  follows: 
s  - 

V'  M  V 

Xs  ”  (V’  M  V>  (A. IV. 3) 


The  position  of  the  average  centroid  x  for  slow  AGC  is  given 

s 

by : 


2. 


Average  Centroid  for  Slow  AGC 

(V’  M  Vx> 

Xs  "  <V"  M  V> 


(A.  IV.  4) 


follows: 


The  operator  (F^ , $2 , . . .  >  is  defined,  as 

2tt  2ir 

<F> '  7777  /■"  /  d*id*2--'<,V<*i-*2 . V 


(A. IV. 5) 


o  o 


The  denominator  of  eqs.  (A. IV. 3)  and  (A. IV. 4)  is  equal  to  the  mean 
radar  cross  section  a  which  is  given  by  eq.  (3a)  of  Section  7.1: 


a  « 


I 


ai 


(3a) 


Combining  eqs.  (A.IV.l)  to  (A. IV. 5)  and  substituting  eq. 


(3a)  gives: 


(A.IV.6) 
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This  formula  is  idsntical  with  eq.  (6)  of  Section  7.1.2. 
All  the  results  of  Section  7.2.2  entitled  "Radar-Target  Glint"  were 
derived  from  it  by  applying  convent ional,  elementary  statistics 
combined  with  Rice's  theory  of  zero-crossings,  as  expounded  in 
Appendix  A. III. 

It  should  be  emphasized  here  that  the  mathematical 
expressions  for  the  effects  of  AGC,  as  given  by  eqs.  (A.IV.31 
and  (A.IV.4),  are  simplified  formulations  of  the  actual  AGC 
response.  A  more  rigorous  formulation  must  account  (1)  for  the 
fact  that  the  AGC  filter  has  a  dynamic  response  given  by  an 
operator  or  a  transfer  function,  (2)  for  the  real  gain  control 
which  is  not  the  Inverse  of  the  AGC  filter  output  but  is 
represented  by  a  more  complicated  function. 


3. 


Instantaneous  Centroid  for  Fast  AGC 
A  simplified  expression  for  the  position  of  the 


instantaneous  centroid  x?  under  the  action  of  fast  AGC  is 
g iven  by : 

V  U  Vx 

xf  “  vnrr 

Substituting  eqs.  (A.1V.1)  and  (A. IV. 2)  in  (A. IV. 7)  gives: 


(A. IV. 7) 


XF* 


fa  ^*1  008  (V*j: 


(A. IV. 8) 


4.  Average  Centroid  for  Fast  AGC 

The  position  of  the  average  or  effective  centtold 
under  the  action  of  fast  AGC  becomes: 


(A.  IV.  9) 
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Combining  eqs.  (A. IV. 5)  and  (A. IV. 8)  one  obtains: 

2w  2v 


„  „  Ll  Ji  xi  Soy^T  cos (#.-♦.) 

— — -  (A-1V- 

kMW 


<2v) 


9a) 


cos  <*i-v 


This  multiple  integral  is  difficult  to  evaluate  and  only 
partial  results  are  known  to  this  author. 

James  Hanson  and  collaborators  of  APL  investigated 
problems  of  this  kind.  In  one  study  they  applied  the  following 
restrictions: 

(a)  The  analysis  is  linearized.  Host  angles  are 
very  small. 

(b)  All  the  are  variable  but  correlated.  They 

ckange  as  linear  functions  of  time  and  stay  in 
commensurate  ratios. 

Then  one  gets  the  following  result: 

X1  1  xp  —  XN  ' 

The  present  author  attacked  this  problem  in  a  different 
way.  Instead  of  using  vectors  and  matrices,  he  cast  the  expression 
(A. IV. 9)  in  the  form  of  complex  variables.  Then  the  following 
result  could  be  easily  proved: 

If  there  exists  a  dominant  component  scatterer  S 

n 

with  overpowering  radar  cross  section  aQ,  i.e.,  for 

°n  »  °1  +  °2  +  •••  +  Vl  +  °n+l  +  •••  +  (A.IV.IO) 


one  gets 

XF  "  xn  ’ 


(A.  IV.  11) 
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This  fact  is  well  known  for  N  -  2.  If  a  radar  applying 
"fast”  or  "instantaneous"  AGC  tracks  two  unequal  targets  within  the 
linear  regime  of  its  bean  width,  then  the  stronger  target  becomes 
dominant  and  pulls  the  radar  axis  to  its  own  position.  For  three 
or  sore  unequal  targets  this  fact  does  not  hold  true,  unfortuuately . 
Inequality  (A. IV. 10)  indicates  how  strong  one  single  target  has  to 
be  in  order  to  pull  the  radar  axis  into  its  own  position. 

If  one  deals  with  a  target  that  consists  of  a  dominant 
scatterer  and  a  series  of  minor  echo  sources,  it  often  seems 
desirable  to  hit  the  dominant  scatterer.  If  this  is  the  preferred 
tactic,  then  the  use  of  fast  AGC  appears  advantageous. 

Muchmore  (Ref.  [l])  investigated  the  effects  of  AGC 
on  radar  target  scintillation.  The  author  assumed  the  following 
conditions: 

(a)  Simplified  target  model  (one-dimensional,  uniform 
echo  density) , 

(b)  Small  error  angles, 

(c)  No  target  maneuvers., 

(d)  No  amplitude  noise,  no  receiver  noise,  no 
jamming  noise. 

Under  these  conditions  Uuchmore  arrived  at  the  conclusion 
that  alow  AGC  is  superior.  He  showed  that  very  fast  AGC  may  increase 
the  scintillation  noise  density  by  a  factor  of  approximately  three. 

Of  course,  if  these  conditions  are  relaxed,  a  different 
situation  will  arise.  In  the  example  given  by  the  present  author 
above,  condition  (a)  is  rescinded.  The  craft  that  possesses  one¬ 
dimensional,  uniform  echo  density,  is  replaced  by  a  target  that 
contains  a  dominant  scatterer  and  minor  echo  sources.  In  this 
case,  fast  AGC  is  advantageous. 

References  [7],  [23]  and  [24]  investigated  the  effects 
of  AGC  under  general  conditions  where  restrictions  (a) ,  (b) ,  (c) , 
and  (d)  are  partially  or  completely  relaxed.  All  three  favor 
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fast  AGC.  Their  conclusion  is  that  under  practical  tracking 
conditions  a  fast-acting  AGC  will  give  better  tracking  performance. 

However,  if  one  can  afford  a  sophisticated,  flexible 
system,  it  might  be  worth-while  to  adapt  the  AGC  to  the  instan¬ 
taneous  situation.  For  instance,  during  the  terminal  flight  of 
a  missile  which  is  homing  on  a  long  target,  glint  might  well  be 
the  dominant  noise  source.  If  so,  it  would  pay  off  to  switch  the 
automatic  gain  control  to  a  setting  with  a  longer  AGC  time  constant. 
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A. V.  ADDITIONAL  NOTES  CONCERNING  OUR  ANALYSIS 

1.  Effects  of  Vibrations 

In  the  derivation  of  the  electrical  phase  angles  of 

the  individual  scatterers  (see  eq.  A. I. 2a  and  A. I. 2b)  we  followed 
the  analyses  of  References  [5]  and  [19].  In  the  monostatic  case 
on'1  gets: 

=  -2k[(Xj-x)  cos  9  -  (y.  -y )  sin  0]  +  $  (A.  1. 2a) 

Some  authors  (see  References  [1],  [21],  and  [22])  modify  the 
expressions  for  the  individual  phase  angles  slightly,  as  follows: 

=  ~2k[(Xj-x)  cos  9  -  (y^-y)  sin  0]  +  a ^  (A.V.l) 

where  the  a ^  are  statistically  independent  random  variables 
that  follow  a  uniform  distribution: 

p (a.)  -  for  0  <  a.  <  2ir  (A.V.2) 

J  eir  —  j  — 

In  the  target  model  by  J.  J.  Freeman  (see  Ref.  [21]) 
the  random  phases  are  supposed  to  account  for  structural 

vibrations  and  possibly  other  random  phenomena.  The  same  author 
made  a  very  diligent  study  of  the  effects  of  these  random  angles 
on  tracking  performance.  One  of  his  results  is  that  their  effect 
on  a  number  of  mean  statistical  output  quantities  (such  as  the 
effective  tai ^et  center,  correlation  functions  and  spectra)  is 
negligible.  This  result  is  gratifying  to  our  studies,  inasmuch 
as  we  did  not  include  the  important  effects  of  structural  vibrations 
in  the  preceding  analysis.  However,  it  can  be  easily  shown  that  if 
we  had  introduced  expression  (A.V.l)  containing  the  random  phases 
ofj  in  our  analysis,  instead  of  equations  (A.  1. 2a)  and  (A.  1. 2b), 
we  would  have  wound  up  with  the  same  results.  For  instance,  in 
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the  example  of  Section  A. II,  taking  the  derivative  of  the  phase 

function  and  equating  it  to  zero  (leading  to  eq.  (A. II. 12))  would 

have  eliminated  the  constant,  random  phases  (a.  -a.  ).  In  the 

ic  jo 

example  of  Section  A. Ill,  selecting  the  subset  (A. III. 4)  of  terms 

and  adding  would  have  canceled  the  terms  of  random  phases 

[(a, -a.)  -  (a  -a) ]  in  eq.  (A. III. 3). 
j  l  m  n 

This  fact  strengthens  our  argument  for  the  application 
of  Kelvin's  principle  and  its  extensions  in  our  analysis,  since 
we  used  an  approach  to  the  problem  which  is  quite  distinct  from 
the  method  of  Reference  [211.  Yet  we  arrived  at  the  same  result 
in  the  case  of  the  effects  of  vibrations. 


2 .  Remarks  Concerning  Soectra 

As  explained  in  Section  II,  spectra  are  not  included  in 
our  fin'll  results  because  they  are  characteristic  properties  of 
stationary  processes.  We  used  the  concepts  of  Rice's  noise  theory, 
such  as  the  average  count  of  zero  crossings  and  other  instantaneous 
statistical  quantities  which  are  applicable  to  locally  stationary 
processes.  Since  we  only  needed  the  ratio  of  spectra  in  the  inter¬ 
mediate  results,  we  could  simplify  our  analysis  by  making  full  use 
of  assumptions  (2)  and  (3)  of  Section  VI.  We  idealized  the  radar 
echo  signals  by  assuming: 


L/N 

X 


CD 


(A.  V.  3) 


where  L  is  the  length  of  the  target,  N  is  the  number  of 
scatterers  and  A.  is  the  RF  wavelength.  As  a  consequence  of 
the  limiting  relation  (A.V.3),  the  spectra  of  the  radar  echo 
signals  became  wide-banded  and  flat, similar  to  the  spectra  of 
"white  Gaussian  noise  processes." 

In  contrast  to  our  wdrk,  References  [1],  [21]  and  [22] 
computed  spectra  as  part  of  their  final  results.  Muchmore 
(Ref.  [1])  applied  the  conditions  enumerated  in  Appendix  A. IV,  pg.  52, 
and  computed  for  the  normalized  spectrum  of  the  IF  target  signal 
the  following  expression: 
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wIF(f) 


1 

2f 


m 


for  f  -f  <  f  <  f  +f 
o  m  —  —  o  m 


0  otherwise 


The  normalized  video  spectrum  becomes: 


WV<1) 


F-  <1  -  2J-)  lor  0  <  f  <  2  im 
m  m 

0  otherwise 


The  scintillation  bandwidth  f  is  given  by: 


f  -  r-  0  cps  . 
m  a 


L  is  the  length  of  a  slender, strip-like  target  which  is 
positioned  at  right  angles  to  the  line-cf -sight  and  turns  with 
a  constant  rate  of  fi  radians  per  sec,  fQ  is  the  carrier 
frequency. 

Freeman  (Ref,  [21])  computed  the  spectrum  of  the 
angular-error  signal  under  similar  conditions  except  that  in 
his  case  the  probability  is  1/2  that  the  target  turns  with  a 

positive  rate  (+n)  and  1/2  that  it  rotates  with  a  negative 

,  * 
rate  ( — Q) .  The  spectrum  in  normalized  form  becomes: 

W£(f)  - 

Again,  the  scintillation  bandwidth  before  the  filtering  action 
of  the  servo,  is  given  by  (A.V.6).  There  is  no  discrepancy 
between  these  results  and  the  findings  of  our  Section  A. III. 


2|~  [1  +  (1-  y~)31  for  0  <  f  <  2  f 
m  m 

0  otherwise 


This  result  is  identical  with  eq.  (21)  of  Muchmore  (Ref.  [1]) 


-  56  - 


(A.V.4) 


(A.V.5) 


(A.V.6) 


(A.V.7) 
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If  tends  toward  infinity  (relationship  A,V.J>,  a  fortiori  In¬ 

tends  toward  infinity.  Hence  the  bandwidths  f  of  ihe  afore- 

m 

mentioned  spectra  due  to  eq.  (A.V.6)  grow  over  all  bounds,  in 
agreement  with  eq.  (A,  III. 5),  computed  for  N  -  oo  . 

It  is  interesting  to  note  that  the  average  time  rate 
of  counts  is  simply  related  to  f^  by: 

N+  ~  c  •  f  (A.V.8) 

t  m 

where  c  is  a  function  of  the  target  configuration,  of  the 
pertinent  variables  (0,  f,  t)  and  of  the  quantity  that  is  counted. 

If  the  target  resembles  the  model  of  references  [1],  [21]  and  [22] 
(one-dimensional,  uniform  distribution  of  scatterers) ,  then  the 
constant  c  can  be  readily  computed.  This  job  is  being  carried 
out  at  the  present  time  and  will  be  documented  in  a  separate 
report . 


3.  The  Use  of  a  Square-Law  Detector 

The  present  analysis  (see  assumption  (8)  of  Section  VI) 
and  most  pertinent  References  ([1],  [61,  [21],  [22])  use  square- 
law  detectors  in  their  investigations  of  radar  target  scintillation. 
Since  many  modern  radar  receivers  apply  linear  envelope  detectors, 
it  is  well  to  dedicate  a  few  thoughts  to  the  implications  of  this 
procedure. 

The  situation  can  best  be  summed  up  by  a  statement  of 
M.I.  Skolnik  (Reference  [26],  pg.  431)  and  we  quote  verbatim: 

"In  general,  the  difference  between  the  two  (detectors)  is  small 
and  the  detector  law  in  any  analysis  is  usually  chosen  for 
mathematical  convenience."  This  fact  is  borne  out  by  the  mathe¬ 
matical  derivation  and  graphical  description  of  the  continuous 
parts  of  the  spectra  of  the  two  detector  outputs,  in  Reference  [25]. 
Comparison  of  Figures  3.13  and  3.14  of  this  reference  indicates 
that  there  is  indeed  little  difference  in  the  shapes  of  the  two 
spectra. 
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Similar  to  our  discussion  regarding  the  selection  of  the 
AGC  (Sec.  A. TV),  the  choice  of  the  type  of  detector  is  dictated  by 
the  tactic j.1  situation  and  the  noise  environment.  If  one  can  afford 
a  sophisticated,  adaptive  system,  a  more  detailed  analysis  is 
warranted.  The  following  remarks  are  added  with  this  objective 
in  mind. 

Delano  (Ref.  [6])  made  a  study  for  the  case  where 
angular  scintillation  is  dominant.  He  obtained  the  following 
result:  Let  T]  be  the  rms  centroid  fluctuation  about  the  mean 
radar  center  for  a  system  using  a  linear  envelope  detector  and  T| * 
the  equivalent  quantity  for  a  square -law  detector.  Then 


(A.V.9) 


This  improvement  of  12%  for  a  square-law  detector  applies  to  the 
final  stages  of  missile  flight  where  target  glint  might  become 
dominant . 

What  is  the  situation  at  the  early  stages  of  flight, 
say  during  detection  and  acquisition?  J.  I.  Marcum  (Ref.  [3], 
pg.  211)  made  a  study  of  the  optimum  detector  law.  For  the 
early  stages  of  flight  where  receiver  noise  is  dominant  and  the 
S/N-ratio  is  small,  he  found  that  a  square -law  detector  closely 
approximates  the  ideal  device.  In  contrast,  for  the  cases  where 
the  S/N-ratio  is  large  (say  during  later  stages  of  tracking)  a 
linear  envelope  detector  gives  ideal  performance.  This  result 
agrees  with  the  findings  of  R.  A.  Smith  (Reference  [28]).  Let 
Pjp  and  Py  be  the  signal -over -noise  ratios  for  the  IF-stage 
and  video-stage  respectively,  using  a  linear  detector,  and  let 
PjF  and  Py  be  the  equivalent  terms  for  a  square -law  detector. 
Then  for  large  signal-over-noise  ratios,  i.e.  ,  pIF  =  pjF  -•  a> 
he  finds: 


2  and 


1 

2  ' 


(A.V.10) 
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The  improvement  by  a  factor  of  4  in  favor  of  the  linear 
detector  is  not  as  significant  as  it  appears,  because  for  noise¬ 
less  signals,  the  type  of  detector  becomes  irrelevant. 

However,  in  a  practical  problem,  the  specific  noise 
environment  should  be  carefully  scrutinized  and  the  system 
designed  accordingly. 


4 .  Remarks  Concerning  the  Counting  of  Lobes 

The  mean  count  of  lobes  or  the  average  rate  of  nulls 
was  defined  in  Section  A. II  for  a  Gaussian,  locally  stationary, 
unbiased  process  by  the  average  number  of  zero  crossings  with 
positive  slope  per  unit  of  the  independent  parameter  and  given  by 
Rice's  formula  (A. III. 6).  If  the  process  is  biased  such  as  the 
uni— polar  video  signal  of  the  target  echo,  the  formula  can  still 
be  used  but  must  be  applied  to  the  unbiased  process  which  is 
formed  by  subtracting  the  mean  value  from  the  biased  process. 

(See  Figure  3) . 


However  in  some  cases  (see  Figure  4),  eq. (A. III. 6)  will 
not  give  the  count  of  lobes  or  nulls  properly.  If  one  is  interested 
in  the  number  of  lobes  or  nulls  in  a  strict  sense,  then  one  should 
count  the  maxima  or  minima  of  the  process  z(q).  Eq.  (A. III. 6) 
has  to  be  replaced  by  another  formula  of  Rice  (cf.  eq.  (3.6-6)  of 
Ref.  [10]).  One  gets  for  the  average  number  of  maxima  per  unit 
of  parameter  q  (say  1  sec  or  1  rad  or  1  Hertz)  the  term  M  (z;q) : 

t  \  *1 


Mq(z;q) 


JL 

2l T 


pz(4)  fr;q) 

-Pz  (t  ;  q) 


1/2 


T-0 


(A.V.ll) 


/ 


For  specific  targets  (A.V.ll)  will  lead  to  expressions  similar  in 
appearance  as  eqs.  (10)  to  (12)  of  Section  VII.  However  the 
mathematical  and  numerical  computations  will  be  slightly  more 
complicated. 
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Mean  Number  of  Zero-Crossings  per  Unit^  of  Parameter  q  of 
Unbiased  Process  Az  =  z  -  z.  Here  N  (Az;q)  =  10  (zeroes 
with  positive  slope) .  q 


Mean  Number  of  Maxima  (Lobes)_  per  Unit  of  q  of  Unbiased 
Process  _Az  -  z  -  z.  Here  N  (Az;q)  =»  6  (zeros  with  positive 
slope)  Mq(Az;q)  -  11  (Maximaqor  Minima). 
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If  the  theory  is  applied  tc  slender  (one-dimensional) 
targets  consisting  of  many,  independent  scatterers  following  a 
random,  uniform  distribution,  eq.  (A.V.ll)  yields  a  quantity 
which  is  roughly  50%  greater  than  the  quantity  computed  by 

eq,  (10)  of  section  VII 

On  the  other  hand,  the  latter  result  was  based  on  a  flat 
video-spectrum  or  a  Dirac  delta-function  correlation  (see  eq.  (A. III. 5) 
of  Appendix  A-III).  References  [l],  [21],  and  [22]  showed  that  the 
video  spectrum  can  be  taken  to  be  triangular  (cf.  eq.  (A.V.5)).  In 
many  cases  it  is  tapering  off  even  more  sharply.  If  our  results 
are  adjusted  for  this  effect,  one  has  to  reduce  M  by  roughly 
33%. 

Combining  these  two  adjustments  for  slender  targets  with 
many  independent,  uniformly  distributed  scatterers  one  obtains  the 
following  result:  Eq.  (10)  of  section  VII  yields  an  N  which  in 
practice  yields  the  count  of  maxima  (lobes).  This  fact  is  born 
out  by  experimental  measurements.  Details  of  analytical,  numerical 
and  experimental  results  will  be  documented  in  future  reports. 
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